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Abstract 

Bose-Einstein correlations in relativistic heavy ion collisions are examined in a general model containing 
the essential features of hydrodynamical, cascade as well as other models commonly employed for describing 
the particle freeze-out. In particular the effects of longitudinal and transverse expansion, emission from 
surfaces moving in time, the thickness of the emitting layer varying from surface to volume emission and 
other effects are studied. Model dependences of freeze-out sizes and times are discussed and compared to 
recent Pb + Pb data at 160A-GeV. 



I. INTRODUCTION 



Bose-Einstein interference of identical particles or the Hanbury-Brown & Twiss effect (HBT) |jj shows up 
in correlation functions of pions and kaons emitted from the collision zone in relativistic heavy ion collisions. 
It is an important tool for determining the source at freeze-out and recent data ^-^| can restrict the rather 
different models, that have been developed to describe particle emission in high energy nuclear collisions. In 
hydrodynamical calculations particles freeze-out at a hypersurface that generally does not move very much 
transversally until the very end of the freeze-out 0. In cascade codes the last interaction points are also 
found to be distributed in transverse direction around a mean value that does not change much with time 
§-jllJ, but the width of the emission zone increases from narrow surface emission to a widespread volume 
emission. Other models like cylindrical symmetric models with Bjorken longitudinal scaling and volume 
emission |l^,[l3| or surface emitting sources 0Jl^] have also been studied and sizes, freeze-out times, etc. 
have been estimated by comparing to experimental data. 

In this letter we want address the general dynamical features of particle emission in relativistic heavy ion 
collisions. We start from a general source, incorporating volume as well as surface emission, longitudinal 
and transverse flow and expansion as well as moving freeze-out surfaces. The various contributions to the 
HBT radius parameters are calculated and we show, that they can be interpreted as fluctuations in radial, 
temporal, angular and emission layer thickness variables. Finally, we discuss a more quantitative analysis 
and compare to recent NA44 HBT data on central Pb + Pb collisions at 160A-GeV. 

For the correlation function analysis of Bose-Einstein interference from a source of size R we consider two 
particles emitted a distance ~ R apart with relative momentum q = (ki — IC2) and average momentum, 
K = (ki + k 2 )/2. Typical heavy ion sources in nuclear collisions are of size R ~ 5 fm, so that interference 
occurs predominantly when q^h/R ~ 40 MeV/c. Since typical particle momenta are ki — K ~ 300 MeV, 
the interfering particles travel almost parallel (see Fig. (1)), i.e., ki ~ fc 2 ~ K ^> q. The correlation function 
due to Bose-Einstein interference of identical particles from an incoherent source is (see, e.g., [p"3"[) 

C2(q ' K) - 1 ± 1 J<fixS(x,K) 1 ' (1) 

where S(x, K) is a function describing the phase space density of the emitting source. The +/— refers to 
boson/fermions respectively. 

Experimentally the correlation functions for identical mesons (n^n^ , K^K , etc.) are often parametrized 
by the gaussian form 

C 2 (q s , q a , m) = 1 + Aexp(-gX 2 - <?X - Q?Rf - ZlodlRli) • (2) 

Here, q = ki — k 2 = {q s ,qo,qi) is the relative momentum between the two particles and Ri,i = s,o, / the 
corresponding sideward, outward and longitudinal HBT radius parameters respectively. We will employ the 
standard geometry where the longitudinal direction is along the beam axis and the outward direction is 
along K and the sideward axis is perpendicular to these. Usually, each pair of mesons are lorenzt boosting 
longitudinal to the system where their longitudinal momentum or rapidity vanish, Y = 0; here their average 
momentum K is perpendicular to the beam axis and is chosen as the outward direction. In this system the 



pair velocity /3jf=K/ 'Ek points in the outward direction with (3 = p±/m± where m± — y/ m 2 + p \ is the 
transverse mass. As pointed out in [jl3| the out-longitudinal coupling R D i only vanish to leading order when 
Y = Q. The reduction factor A in Eq. (|^) may be due to long lived resonances ]l8| , p^[ , coherence effects, 
incorrect Gamov corrections or other effects. It is found to be A ~ 0.5 for pions and A ~ 0.9 for kaons. 




FIG. 1. Cross section of the interaction region perpendicular to the longitudinal or z-direction. Particles have to 
penetrate a distance I out to the surface of the interaction region in order to escape and reach the detector. The bulk 
part of the emitted particles comes from a surface region of width X m fp, the mean free path of the particle. 

It is very convenient to introduce the the source average and fluctuation or variance of a quantity O defined 

by 

With qx ~ q • x — q-/3x t one can, by expanding to second order in Ri and compare to Eq. (g), find the 
HBT radius parameters Ri, i=s,o,l. ^ They are (see, e.g. |l3|]) 

Ri = <r(z 4 -Pit). (4) 



"The expansion at small qi and extraction of the HBT radius parameters Ri, Eq. (^), may not be directly comparable 
to the gaussian radii extracted experimentally because the experimental data has best statistics around qi ~ h/Ri ~ 
40MeV/c. For the specific models discussed here we have checked that they are comparable. We attribute this 
agreement to the many different fluctuations contributing to the HBT radius parameters which "smear" out the 
source; for many random fluctuations we expect that the centrahlimit theorem guaranties a gaussian distribution. 
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The HBT radius parameters are a measure for the fluctuations of (xi — Pit) over the source emission function 
S. 

Particle production in ultrarelativistic heavy ion collisions has been found to have dynamical features of 
strong longitudinal expansion and some transverse expansion p0| . Furthermore, the geometry plays a role 
and one would expect that particles escape from the outer layers of the fireball or if the source is sufficiently 
small the whole volume will freeze-out. We model such a general cylindrically symmetric source by assuming 
local thermal equilibrium with longitudinal Bjorken flow (u z = z/t) as well as transverse flow v through a 
Boltzmann factor. The space-time geometry is modeled through the transverse source distribution S±(r±_) 
and the temporal part S t (t) 

S{x,K) ~ e- K - u ' T S ± (r ± )S T (r). (5) 



Here, r = \/t 2 — z 2 is the invariant time and r/ = 0.51n(t + z)/(t — z) the space-time rapidity. Including 
transverse flow v = (v x ,Vy) the flow four-vector is u — 7(w)(cosh(ry), sinh(ry), v) [r^Jl3] |, which gives 



K ■ u — m^7(u)(cosh(?7 — Y) — O ■ v) . (6) 

Here m± is the transverse mass and Y the rapidity of the particles. In the following all variables are boosted 
into the frame in which Y = 0. Notice that any normalization cancels out in correlation function (0). One 
can in addition apply a factor S v (rf) ~ exp(— (rj — r] cms ) 2 /25r] 2 ) to correct for the lack of Bjorken scaling 
near target and projectile rapidities [^2[ . However, the thermal factor exp(— m± cosh(?7 — Y)/T) centers 
the space-time rapidity rj around the pair rapidity Y on a scale ~ y/T '/m±_ much narrower than Srj. The 
additional factor has thus only a minor effect, reducing the longitudinal source sizes slightly. 

It is important to distinguish between volume and surface freeze-out, since they give very different HBT 
radius parameters [lq| . Hydrodynamical models |?J assume that particles are emitted from the surface. This 
is actually also found in some cascade models at early times of the collision |8|-[ll| , but eventually the whole 
source freezes out and disintegrates. The late stage of cascade models resembles more a volume freeze-out. 
We want to mimic these different models and stages by introducing a source which emits particles according 
to a simple Glauber theory. The surface freeze-out component is described by a Glauber absorption factor, 
which suppresses particles escaping from the interior of the source. The source becomes opaque 

Sl (rj- e -S? dx '° p{x,) , (7) 

where the integral runs over the particle trajectory from last interaction point x. Modifications of single 
particle spectra in hydrodynamic calculations due to such an emission layer (0) has been consider in Jl^ ] . We 
introduce the mean free path X m f p = (cr(p)) -1 , where (p) is the average density in the emission layer and a 
the interaction cross section. Glauber absorption only allows emission from a layer of thickness ~ X m fp just 
inside the surface radius R. Thus we can rewrite 

S±(r ± ) ~ e' [ / x ^Q(R-r ± ), (8) 



where I = \/ R 2 — y 2 — x is the distance the particle has to pass through the source in order to escape from 
the surface when (x, y) = R (cos#, sin#) is its position in outward and sideward direction respectively (see 
Fig. (1)). In hydrodynamical calculations \ m fp = but in cascade codes it can be several fm's . The surface 
may also move in the transverse direction with time, R(t). 

The temporal emission of the source is determined by S T (r) . It is commonly approximated by a gaussian, 
S t (t) ~ exp(— (r — tq) 2 /2(5t 2 ), around the source mean life-time, tq with width, St, which is the duration 
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of emission. These gaussian parameters approximate the average emission time, (r) and the variance or 
fluctuation, cr(r), respectively for a general source. 

An opaque source emits from a thin surface layer instead from the whole volume. We can thus calculate the 
HBT radius parameters very generally by expanding in X m fp/R 1- For strict surface emission, X m f p — 
and the source (||) reduces to 

Sj_ ~ 8{R{t) -r_L)6(cos0) cos0. (9) 

The geometric factor cos 9 suppress the peripheral zones and the 0(cos#) factor insures that particles are 
only emitted away from the surface, i.e., only particles from the surface layer of the half hemisphere directed 
towards the detector will reach it whereas particles from the other hemisphere will interact on their passage 
through the source. t Most hydrodynamical freeze-out mechanisms do not include the directional condition 
that particles can only be emitted away from the surface though strong flow at the surface has a similar 
effect. Whereas the source in @ corresponds to black body emission, which is constant in time per surface 
element, the temporal variation is described by S t (t). 

The HBT radius parameters can now be calculated from Eq. (|J) using Eqs. (JsJ> and (||) for any mean free 
path and transverse flow and results are shown in Fig. 2. It is, however, very instructive to consider the 
case of emission from a surface layer, X m f p <C R(t), where a number of simplifications appear and analytical 
results can be obtained. Since the emission points are narrowly confined within a layer of thickness A m / P 
within the surface, an expansion in the mean free path can be performed and the integration in transverse 
radial direction simplifies. The space-time rapidity, angular and temporal integrations separates and due to 
the normalization a number of factors cancel when evaluating the HBT radius parameters. For example, a 
function of proper time only needs to be averaged with respect to the temporal parts of the source 

J T UrrR(r)S T (r)O(r) 

U vr)) = Tt]-, , w : — ■ (10) 

Jo tR(t)S t (t) 

The angular averages also simplify for the cylindrical geometry. From the definitions in Eq. (|^) we obtain 

r /2 /0 O(0) exph{v s )v sP± cosO/T) cos 9 d9 

{O(0)) = — — — . (11) 

J_4/ 2 exp(7(u s )w s px cos9/T) cos9d9 

The factor exp(-f(v s )v s p± cos 9/T) includes the effect of transverse flow v s at the surface (see Eq. (||) and 
has the effect of narrowing the angular emission in the direction of K, i.e. in the outward direction (9 = 0). 
A spherical source would have an additional factor R(t) in ( |To| ) and sin0 in the integrals of Eq. (fil|). 

The HBT radius parameters can to leading order in fluctuations and to second order in X m fp *C R(t) (in 
the frame where Y — 0) now be evaluated * 



1 

6' 

?2 _ 



R 2 s = <Y) = (y 2 > = (R(r) 2 )a(sm9) - -X 2 mfp , (12) 



R 2 = (T(X - p t) 

= (R(t)) 2 a (cos 9) + a (R(t)) (cos 2 9) + 0^(7 (t) 



^Detectors on the other side of the beam line measure particle from the other hemisphere. Thus relativistic heavy 
ion collisions have an advantage to stellar interferometry which cannot measure the back side or, as referred to in the 
case of the moon, the "dark side" . 
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-2/3 (co S 6)((R(t)-(R(t)))(t-(t))) + (---)X 2 mfp , (13) 

(t 2 ) T 

R? = a(z — Pit) — f7(rsinh?7) ~ — — — . (14) 

j(v s ) mj_ 

The terms in R 2 are the average of the square of the transverse source size times the angular fluctuations 
which, due to the inversion symmetry along the y-axis, is a(sm9) = (sin 2 9)^ As the emission layer is inside 
the source particle have r± < R and R s the sideward HBT radius parameter is reduced by a term of order* 

\2 
rnfp- 

The terms in R 2 are respectively: i) angular fluctuations* : ct(cos6>) = (cos 2 9) — (cos#) 2 , ii) fluctuations in 
transverse radial direction: <t(R(t)) = (R 2 (t)) — (i?(r)) 2 , in) temporal fluctuations: <j(t) = (r 2 ) — (r) 2 , iv) 
a cross term between radial and temporal variations which is positive for an inward moving surface, v) and 
finally a term of order X 2 n ^ p due to the thickness of the surface layer adding to the outward fluctuations. In 
addition the average life-times and fluctuations in the life-times of short lived resonances should be added 
to the temporal fluctuations |lS[j ; the long lived resonances can account for most of the reduction factor A 
in Eq. (§). 

It was assumed that Y = 0, i.e. the longitudinal pair velocity vanishes (3i = 0. The longitudinal HBT 
radius parameter is then the fluctuation o~(z) — cr(rsinh?7). When T <C m± the space-time rapidity n is 
small and one finds cr(sinh?7) ~ o~(rj) ~ T/m±. 

A transparent source corresponds to the opposite limit of (||), i.e. A m / P — > oo. The transverse radius 
parameters are in this case the same R s = R = R/2 for a source without transverse flow or fluctuations 
in (duration of) emission time. The longitudinal HBT radius parameter is unchanged for m± 3> T, since 
it is only sensitive to the longitudinal expansion and emission time but not to transverse coordinates. The 
general dependence on transparency/opacity is shown in Fig. (2), where the transverse radius parameters 
are plotted as functions of the transverse flow parameter ~/v s pj_/T for different X m f p /R. A flow profile of 
the form v = v s t/Rq was assumed. For simplicity we have neglected the radial dependence hidden in "f(v s ), 
which is allowed for small mean free paths or small transverse flow velocities. For transparent sources it 
reduces the transverse HBT radius parameters by the same amount to order v 2 Jl^| . Also notice that only 
spatial fluctuations are included in Rq whereas temporal fluctuations such as (3 2 <r(r) still have to be added. 



*For a cylindrically symmetric source without transverse flow o(sm6) — (sin 2 8} — 1/3, {cos 2 8} — 2/3, and 
cr(cos<9) = 2/3 — (7r/4) 2 — 0.05. In the coefficients to A^ p we have, for simplicity, neglected transverse flow. 
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For a transparent source with A m / P S> R and without transverse flow the outward HBT radius parameter 
is larger than the sideward p2| , p^ | 

Rl = R 2 s +(3 2 <j(r). (15) 

The excess is due to the duration of emission, cj(t) 1 / 2 , of the source in which particles with outward ve- 
locity p o = p±/mj_ travel on average a distance /JqO^t) 1 / 2 towards the detector. The sideward distance is 
perpendicular to this velocity and R s is therefore affected by the duration of emission and thus reflects the 
"true" transverse size of the source. The longitudinal HBT radius parameter Ri is affected by the duration 
of emission through the factor (t 2 ) = (t) 2 + <j(t). 

In relativistic heavy ion collisions the outward and sideward HBT radius parameters are measured to be 
similar [^]-|| and in a few cases the outward is even measured to be smaller than the sideward HBT radius 
parameter ^,|| contradicting Eq. (p"5|), however, within experimental uncertainty. According to Eq. ( |l5| ) 
this implies that particles freeze-out suddenly, St -C Ri, as in a "flash" @, in particular when resonance 
life-times are included Jl8|,[l2|]. However, both the opacity effect leading to surface emission and transverse 
flow reduces R more than R s and so it is possible that R < R s . The opacity effect is evident from Eqs. 
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( |l2| ) and ( |l3| ) and is due to the surface emission which narrows the emission to (x 2 ) ~ (x) 2 (see Fig. (1)) 
which reduces the fluctuations c(x) significantly. Transverse flow narrows the emission along K, i.e. in the 
outward direction, as can be seen from Eq. (11) and Fig. 2 which therefore reduces R s more than R Q . For a 
source with surface emission and/or transverse flow it is therefore possible to obtain R a < R s but a number 
of other effects will add fluctuations such as <t(t), <t(R(t), \ 2 n f p and more will add on i? 2 (see Eq. (p"7|). 

To investigate the possibility R a < R s it is instructive to consider a source similar to what is found in 
hydrodynamic models. The freeze-out surface remains at an almost constant transverse distance R initially 



but eventually moves rapidly inwards at final time r/. More specifically we choose R(t) = TZ^J 1 — t 2 /t 2 
and constant particle emission per surface element (S t (t) = constant). Inserting the source in Eqs. (|l2|jl4|) 
we find without transverse flow (v = 0) 



R 2 S - -Tl 2 - ^lifp i ( 16 ) 



Rl = Cl K 2 + Cl /3 2 rj + c 2 (3 T f TZ + (I - ^)X 2 mfp ■ (17) 

Rf= 2 -— r|. (18) 
5 mi 

where c x = 2/5 - (3tt/16) 2 ~ 0.053 and c 2 = 2(3tt/16) 2 - vr/5 ~ 0.066. The terms in R 2 a correspond to the 
terms in Eq. (|l^) from angular and radial fluctuations, the temporal fluctuations, the cross term, and finally 
the fluctuations from the width of the surface layer respectively. 

In the recent NA44 data on central Pb + Pb collisions at 160 A-GeV pion HBT radius parameters of 
R s ~ R a ~ 4.5 — 5.0 fm and Ri ~ 5 — 6 fm are found The average transverse momentum was p± ~ 165 
MeV such that O ~ 0.76. ^From the transverse momentum slopes of pions, kaon, protons and deuterium 
in |2FJ] one finds a temperature T ~ 120 MeV and transverse flow (v 2 )*- 1 / 2 -* ~ 0.6c. The transverse flow 
effects on the HBT radius parameters (see equation ( |Tl| ) and Fig. 2 with jv s p±/T ~ 1) are small for these 
numbers which allows us to use the HBT radius parameters determined in equations ( ^6[]T^ ). As seen from 
equation (|l^) and (|l7]) and Fig. 2 the sideward radius is relatively less affected by the Aj^y correction than 
the outward. We extract a freeze-out time rt ~ 12 fm/c and an initial transverse source size 1Z ~ 11 fm from 
the experimental Ri and R s respectively. From the experimental value for R Q we finally extract A m / P ~ 2 — 3 
fm. The initial transverse source size is larger than the geometrical size of Pb ~ 7 fm by an amount in excess 
of the uncertainty in the impact parameter, expected from the centrality cuts. Some expansion of the source 
seems to have taken place before final freeze-out. The average emission time is J (r 2 ) = y / 2/5r/ ~ 7 fm/c 
which is entirely consistent with the freeze-out time required to explain the enhancement in the it~ /tt + ratio 
at low p± due to Coulomb repulsion in the same Pb + Pb collisions [|T| . 

In the above example the pions arc emitted during the whole period from collision to freeze-out and do 
not appear as in a "flash" . The opacity effect is more important in reducing R Q than transverse flow at 
this p±. Fluctuations from radial, angular, temporal, cross term and thickness of emission layer contribute 
by similar amounts to R Q . At large transverse momenta the sideward and outward HBT radius parameters 
are reduced by transverse flow and the longitudinal radius parameter scales like 1 /mj_ . These results are in 
qualitative agreement with the experiment Q . 

Eqs. (p^-pTf) leads to R < R s for small f3 and X m f p . It is, however, crucial to realize that this simple 
model breaks down when the inward moving surface speed exceeds the particle velocity outwards [p] . When 
that occurs the particles are overtaken by the surface and experience a volume freeze-out rather than the 
Glauber picture employed above in which the particles scatter their way out through the surface. Near 
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freeze-out the surface moves rapidly inwards and Eqs. ( |i6f|l8| ) are therefore not valid - especially at low 
A> = P±/mj_. 

As mentioned cascade codes results in sources that can approximately be described by two components, 
initially surface emission but eventually volume freeze-out. Generally, for a two-component source 

S(x)= P S 1 (x) + (l-p)S 2 (x), (19) 

with normalization J d A xSi{x) = 1 such that p is the fraction of particles from source 1, the fluctuations in 
a quantity O is by inserting ([lj]) in (||) 

a(0) = P a 1 (0) + (1 -p) (7 2 (0) +p (1 - p) «0)i - {0) 2 f . (20) 

Here, = J d 4 xSi(x)0 and cr^O) = (O 2 ); - (O) 2 . The fluctuations are the weighted sum of the 

fluctuations in the individual sources and an additional cross term. Since (y) = and (z — Pit) also vanishes 
for Y — this additional cross term does not contribute to the sideward and longitudinal HBT radius 
parameters. It is, however, nonnegligible for the outward HBT radius parameter. 

The sideward and longitudinal HBT radius parameters are less sensitive to the details of the freeze-out 
than the outward HBT radius parameter, which allows us to estimate the transverse source size and freeze- 
out. For example, for a source emitting from a static surface we find R s = TZ/VS @, for volume emission 
of constant density we obtain R s = 1Z/2 (see above) and for an inward moving source R s is given by Eq. 
©. These expressions give a transverse source size in the range 1Z *~ 9 — 11 frn from the NA44 data on 
R s . The longitudinal HBT radius parameter depends on y (r 2 ) which is simply Tf for a source emitting 
particles at freeze-out but smaller by a factor y/2/5 for constant emission per surface element, Eq. (|l8|). The 
freeze-out time thus lies in the range r/ ~ 8 — 11 frn/c using the NA44 data on Ri. In contrast the outward 
HBT radius parameter is very model dependent as it contains contributions from a number of unknown 
quantities such as radial, angular, temporal and width of the emission layer, cross terms between radial and 
temporal correlations as seen from Eq. (|l7]) as well as transverse flow. In addition to the strong dependence 
on whether the source is opaque or transparent, a mixture of these sources produces a cross term in i? 2 , Eq. 
(p0|), besides the individual fluctuations. 

In summary we have expressed the HBT radius parameters as averages and fluctuations in spatial and 
temporal quantities for a general class of sources with longitudinal and transverse expansion as well as surface 
and volume emission changing with time. The outward HBT radius parameter R a consists of fluctuations 
in radial, angular and temporal quantities, the width of emission layer as well as life-times of short lived 
resonances. Glauber absorption leads to emission from a surface layer away from the source and has the effect 
of reducing R a significantly but increasing the sideward HBT radius parameter, R s . Strong transverse flow 
reduces both transverse HBT radius parameters at large transverse momentum - in particular R a . Finally, 
we used the recent central NA44 Pb + Pb data to extract HBT radius parameters and to investigate the 
model dependences of the transverse sources sizes. The sideward and longitudinal HBT radius parameters 
are less sensitive to the details of the freeze-out which allows us to estimate the initial transverse source size 
i?o ~ 9 — 11 fm and freeze-out time Tf ~ 8 — 11 frn/c. The outward HBT radius parameter is, however, very 
sensitive to a number of model dependent quantities leading to the above mentioned fluctuations. 
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